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Abstract 



We present the first study of the light hadron spectrum and decay constants for 
quenched QCD using an 0(a)-improved nearest-neighbour Wilson fermion action at 
P = 6.2. We compare the results with those obtained using the standard Wilson 
fermion action, on the same set of 18 gauge field configurations of a 24^ x 48 lattice. For 
pseudoscalar meson masses in the range 330-800 MeV, we find no significant difference 
between the results for the two actions. The scales obtained from the string tension 
and mesonic sector are consistent, but differ from that derived from baryon masses. 
The ratio of the pseudoscalar decay constant to the vector meson mass is roughly 
independent of quark mass as observed experimentally, and in approximate agreement 
with the measured value. 



Introduction The Wilson formulation of lattice fermions leads to matrix elements which 
differ from their continuum counterparts by terms of order the lattice spacing, a. These 
corrections, which cannot be calculated reliably in perturbation theory, represent a substan- 
tial uncertainty in the determination of hadronic matrix elements. A tree-level improved 
theory has been proposed 0], which in perturbation theory eliminates all terms of order 
((^o)"^^ log" a p. Our objective, here, is to look for evidence of such improvement in masses 
and decay constants of the light hadrons. In this letter, we summarise the results of a study 
in quenched QCD on a set of 18 configurations of a 24^ x 48 lattice at /5 = 6.2, using both 
the standard Wilson fermion action: 

S"^ = a^^-|g(a;)g(x) + K^[g(x)(7/,-r)f/^(x)g(x + /i)-g(a: + /i)(7^ + r)f/^(a:)g(a;)] I (1) 
X { ^J. ) 

and the nearest-neighbour 0(a)-improved, or "clover", fermion action 

Sp = - irgoK^"^ ^ q{x)F^^{x)(T^^q{x). (2) 

F^u is a lattice definition of the field strength tensor, which we take to be the sum of the 
four untraced plaquettes in the /xz/ plane with a corner at the point x 
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Full details of the simulations will be presented elsewhere 



(3) 



Computational Details The gauge field configurations and quark propagators were ob- 
tained using the 64 i860-node Meiko Computing Surface at Edinburgh. The full machine 
has 1 GByte of internal memory and, for these programs, sustains around 1.5 Gfiops. 

We update the SU(3) gauge fields with a cycle of 1 three-subgroup Cabibbo-Marinari heat- 
bath sweep followed by 5 over-relaxed sweeps. The results presented here are based on an 
analysis for each action of the same set of 18 configurations, starting at configuration 16800 
and separated by 2400 sweeps. Propagators are calculated for r = 1 at k = 0.1510, 0.1520, 
0.1523, 0.1526 and 0.1529 for the Wilson action, and at k = 0.14144, 0.14226, 0.14244, 
0.14262 and 0.14280 for the clover action; the latter values were chosen to match roughly 
the pion masses computed in the Wilson case. We use an over-relaxed minimal residual 
algorithm with red-black preconditioning for propagator calculations, with point sources 
and sinks. 



String Tension It is useful to have an estimate of the scale from a quantity which is 
independent of the fermion action. For this reason, we present results for the string tension 
obtained from generalised RxT Wilson loops, using a variational scheme comprising 40 and 



28 recursive blockings |Q for the spatial paths. We obtain the potential by extrapolating the 
3 : 2 and 4 : 3 T-ratio results using an estimate of the excited state contribution. 

We find an acceptable fit to the static potential for 2 < R/a < 12 with the conventional 
lattice Coulomb plus linear term: 

V{R)=C-^ + KR (4) 

and obtain E = 0.274(6) and Ka'^ = 0.0259(9), where the quoted errors are purely statistical. 
Using this value for the string tension to set the scale, by requiring ^/K = 0.44 GeV, yields 
= 2.73(5) GeV. These results are in agreement with a similar analysis on 20^ lattices at 
the same /?- value [|, |^. We find from a comparison with lower (3 results, that asymptotic 
scaling of the string tension does not hold at /3 = 6.2. 

Hadron Spectrum We identify the timeslice interval 12 to 16 as a common plateau re- 
gion in the effective mass plots. We obtain the amplitudes and masses of mesons (baryons) 
by correlated least-x^ fits of a single cosh (exponential) function simultaneously to the ap- 
propriate forward and backward propagators. The fits all have x^/dof in the range 0.3 to 
3. The quoted errors are 68% confidence level bounds estimated by a bootstrap procedure. 
The details of our statistical analysis will be presented elsewhere 

The pion, p, nucleon and A mass estimates in lattice units for the two actions are given 
in table Edinburgh plots for the two actions are given in fig. |T[ The plots are broadly 
consistent, showing a trend towards the physical value for niN/mp with decreasing pion mass. 

We now consider chiral extrapolation of the hadron masses in order to compute the lattice 
scale for the two actions. We perform correlated linear fits to the data at all five k values. 
From the bootstrap analysis, we find that there are strong correlations between the pion 
masses at different n values. The behaviour of as a function of 1/2k for both actions is 
completely consistent with PCAC behaviour throughout the quark mass range used. From 
the chiral extrapolation, we obtain 

fi;c = 0.15329+ J (Wilson), = 0.14314+ \ (clover). (5) 

The lattice scales obtained from correlated linear extrapolations of the other hadron masses 
to the chiral limit are presented in table 0. For both actions, the scales derived from rup 
agree well with the scale from the string tension, whereas the scales from the baryon masses, 
whilst consistent with one another, are too low, in common with previous studies using the 
Wilson action. This discrepancy between the baryon and string tension scales is clearly seen 
in the estimate of the probability distributions for the scales obtained from the bootstrap 
samples, shown in figure ||. 
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Table 1: Hadron masses in lattice units. The last row for each action contains the values 
obtained by linear extrapolation to am.,^ — 0. 
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Table 2: Scales determined from different physical quantities. 
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Figure 1: Edinburgh plots. 
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Figure 2: Distributions of the scales determined from nip and mjv for 1000 bootstrap samples 
compared with the scale from the string tension. 



In an attempt to highlight any differences arising from use of the clover action, we present 
estimates for the vector-pseudoscalar meson mass splitting, which in QCD-inspired quark 
models arises from spin interactions and therefore may be corrected at 0{a) by the spin term 
in the clover action. Experimentally, for both light-light and heavy-light mesons, my — im?p 
is very nearly independent of quark mass: 0.57 GeV^ (p — tt), 0.55 GeV^ {K* — K, D* — D, 
B* — B). Taking the scale from the string tension, the range of pion masses for which we have 
data is 330 MeV to 800 MeV. Thus, it is of interest to ask whether our data for m^ — m^, with 
equal-mass quarks, has similar behaviour. In figure |^, we plot the quantity {anipY — {am-i^Y, 
calculated directly from the bootstrap masses, versus {amT^Y for both actions. The Wilson 
data is consistent with previous work |^ which indicated a negative slope, inconsistent with 
experiment at large quark mass (e.g. for the charmonium system — m?p = 0.72 GeV^). 
Apart from a constant vertical shift, the best estimates from our clover data at the four 
highest quark masses agree with the Wilson points; the larger statistical errors in the clover 
data leave open the possibility of a reduced, or zero dependence on quark mass. Using the 
string tension scale, the experimental value corresponds to 0.075 in lattice units, consistent 
with the three lightest Wilson points and all the clover points. 

Meson Decay Constants The pion decay constant is defined through the matrix element 
of the axial current: 



and our normalisation is such that the physical value is 132 MeV. In table ^, we present the 
values obtained for and for the dimensionless ratio /^/m^, using the Wilson action with 
the local axial current, and the clover action with the "improved" axial current: 



Our lattice results for fn/^p vary only slowly with quark mass, in agreement with the 
experimental observation that fn/iTip (0.17) is approximately the same as fx/^K* (0.18). 

In order to determine the physical values, the lattice results given in table |^ need to be 
multiplied by the appropriate renormalisation constant, or Z^. ~ 1 — 0.132(7^ in 
perturbation theory. If we use the bare coupling constant as the expansion parameter, then 
Z'^ ~ 0.87 and we find, for the Wilson action, fn/'mp — 0.18 ^ |. If instead we use the 
"effective coupling" proposed in reference g'^g ^ 1.75(7q, then Z^ ~ 0.78 and we obtain 
fw/fTT'p — 0.16 ^ |. The perturbative estimate of the renormalisation constant Z^ is close 
to 1 {Z^ ~ 1 — 0.02(7^); using the bare coupling leads to Z^ ~ 0.98, whereas the effective 
coupling gives Z^ ^ 0.97, yielding, for the clover action, f-w/rnp ~ 0.12 + 2 both cases. 
We note that the uncertainty due to the choice of the perturbative expansion parameter is 
about 10% with the Wilson action and only about 1% with the clover action. Folding this 



(0|g(0)7M75g(0)|7r(p)) = ^p^ 



(6) 



TQj ■* — TQj 

g(x)(l + y7- Z})7m75(1 - y7- D)q{x). 
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Figure 3: a^m^ — a^m^ versus a^m^, in lattice units; the left- most point in each plot is 
obtained from the chiral extrapolation of the individual masses. 
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Table 3: Values of the pion decay constant, in lattice units, and the ratio fT^/nip. The last 
row contains values obtained by a linear extrapolation to the chiral limit. 
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Table 4: Values of the p decay constant. The last row contains values obtained by a linear 
extrapolation to the chiral limit. 



in, the overall uncertainty in the determination of fn/iTip may be gauged from the bootstrap 
distributions in figure ^ It is apparent that, within errors, the two actions give consistent 
results. 

In table H, we present values of the p decay constant, defined by the relation 

2 

(0|g(0)7^g(0)|p) = -^e^ (8) 
Jp 

where is the p polarisation vector. For the Wilson action we use the local current and 
for the clover action we use the improved local current. Here, the renormalisation constant 
for the Wilson action is given in perturbation theory by Zy" ~ 1 — O.lYAg"^. Using the bare 
(effective) coupling constant, ^ 0.83 (0.71), which leads to 1/fp ~ 0.39 + j (0.33 + \). 
For the clover action, Zy ~ 0.90 (0.83) when the bare (effective) coupling constant 
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Figure 4: Distributions of fT^/nip corresponding to the two choices for and Z'^ described 
in the text, for 1000 bootstrap samples. 



is used, which gives l/fp ~ 0.34 + 2 (0.31 + 2). The physical value of l/fp is 0.28(1). 
The uncertainty in the renormalisation constants makes it difficult to draw &m conclusions 
about improvement from l/fp computed with local currents. The perturbative uncertainty 
would be removed entirely by use of the conserved vector current. However, since the matrix 
element in equation (^ is non-forward, the use of the conserved current would not eliminate 
the 0{a) corrections in the Wilson case. 



Conclusions We have presented the ffist study of the light hadron spectrum and decay 
constants for quenched QCD using an 0(a)-improved nearest-neighbour Wilson fermion 
action aX (3 = 6.2. Having compared the results with those obtained using the standard 
Wilson fermion action, on the same set of 18 gauge field configurations of a 24^ x 48 lattice, 
we see no statistically significant difference between the best estimates for quantities derived 
from 2-point functions, for pseudoscalar meson masses in the range 330-800 MeV. The clover 
data is typically noisier than the Wilson data. Nevertheless, there may yet be advantages 
from using the clover action for 3-point functions, where improved chiral behaviour for the 



B parameter of K decay has been noted |]I0|. We find that aX (3 = 6.2 the scales obtained 
from the meson sector are consistent with that from the string tension, but there is clear 
evidence of an inconsistent scale from the baryon sector, fp/my is roughly independent of 
quark mass as observed experimentally, and the numerical value is broadly consistent with 
the measured value. 
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